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Abstract 

We report the observation of a new dwarf nova, OT J055717+683226, during its first-ever recorded 
superoutburst in December 2006. Our observation shows that this object is an SU UMa-type dwarf nova 
having a very short superhump period of 76.67 ± 0.03 min (0.05324 ± 0.00002 d). The next superoutburst 
| was observed in March 2008. The recurrence time of superoutbursts (supercycle) is, hence, estimated to 

£C) . be ~ 480 d. The supercycle is much shorter than those of WZ Sge-type dwarf novae having supercycles 

of ^ 10 yr, which are a major population of dwarf novae in the shortest orbital period regime (^ 85 min). 
^v^j . Using a hierarchical cluster analysis, we identified seven groups of dwarf novae in the shortest orbital period 

regime. We identified a small group of objects that have short supercycles, small outburst amplitudes, and 
large superhump period excesses, compared with those of WZ Sge stars. OT J055717+683226 probably 
belongs to this group. 

Key words: accretion, accretion disks — stars: novae, cataclysmic variables — stars: individ- 
ual^ J055717+683226) 

5~ i 1- Introduction outburst behavior of SU UMa stars can be explained by 

the thermal-tidal instability model for the accretion disk 
Dwarf novae (DNe) are a subclass of cataclysmic vari- around the white dwarf (for a review, see Osaki 1996). 
able stars (CVs) and consist of a close binary system con- It is widely believed that the evolution of a CV hav- 
taining a white dwarf and a Roche-lobe filling red-dwarf ing P or b ^ 3 hr is driven by angular momentum removal 
star (Warner 1995a). The orbital periods (P rb) of or- from the binary, associated with gravitational radiation 
dinary hydrogen-rich DNe range from 76 min to 5.7 d (Paczynski 1981; Rappaport et al. 1982). By losing an- 
(Ritter, Kolb 2003). DNe having short orbital periods of gular momentum, the CV evolves toward a shorter P or b 
76 min ^P or b ^3 hr tend to exhibit two types of outbursts, regime, until the onset of degeneracy in the secondary 
normal outbursts and superoutbursts. Such systems are star. After the onset of degeneracy, the CV evolves to- 
called SU UMa-type DNe. SU UMa stars exhibit short- ward a longer P or b regime since the mass-radius rela- 
term periodic modulations, so-called supcrhumps, during tion of the secondary star changes. This scenario gives a 
superoutbursts. Supcrhumps can be a good indicator of qualitative explanation for the presence of a short-period 
the orbital period of a binary system because it is well cut-off at P or b ~ 76 min in the observed P or b distribution, 
known that their periods arc slightly (about a few percent) However, the predicted period minimum (P min ) is signif- 
longer than the orbital period (e.g. Warner 1995b). The icantly shorter than the observed minimum (Kolb 1993; 
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Kolb, Baraffe 1999; Renvoize et al. 2002). 

There are two notable exceptions of hydrogen-rich CVs 
having P orb much shorter than the well-known P m i n at 
- 76 min: V485 Cen (P orb = 59.0) and EI Psc (P orb = 
64.2). A characteristic feature of these two is clear TiO 
absorption bands observed in their quiescent spectrum, 
which indicate luminous secondary stars. It is proposed 
that these systems have evolved secondaries and as a result 
their evolutionary paths are different from that of CVs 
having a main-sequence secondary (Augusteijn et al. 1996; 
Thorstensen et al. 2002a). Uemura et al. (2002) proposed 
that V485 Ccn-like objects may be progenitors of helium 
CVs, so-called AM CVn stars, which are recognized as 
being interacting double white dwarfs. If this is the case, 
we can expect that progenitors of V485 Cen-like objects 
are hidden in the CV population having P or b > 76 min. 
Such objects should have evolved secondaries, and hence, 
high mass-transfer rates compared with CVs having main- 
sequence secondaries (Podsiadlowski ct al. 2003). 

Thorstensen et al. (2002b) pointed out that there is a 
large diversity in DNe having very short P or b , particularly 
in terms of the recurrence time of superoutbursts (super- 
cycle; T s ). It is well known that most systems have quite 
long supercycles of T s > 10 yr in the shortest period regime 
of 76 min < P or b ^ 85 min. These are called WZ Sge stars 
(Kato et al. 2001b). The shortest P or b regime also includes 
ER UMa stars, such as DI UMa (P orb = 78.5722 min), 
which have very short T s (20-50 d) (e.g. Kato, Kunjaya 
1995; Ishioka et al. 2001a). V844 Her is a unique object 
also having quite short P or b (78.6859 min) in terms of the 
intermediate T s (260 d), which is rather typical for a DN 
having a longer period of P or b > 90 min (Kato et al. 2003a; 
Kato, Uemura 2000). The nature of the diversity in T s is 
poorly understood. 

According to the disk instability model, T s depends on 
the mass-transfer rate; a higher mass-transfer rate yields 
a short T s (Osaki 1995b). The mass-transfer rate at each 
P or b depends on the binary parameters, such as the struc- 
ture of the secondary star. Hence, the diversity in T s 
possibly indicates that DNe in the shortest P or b regime 
includes several groups with different secondary star struc- 
tures, or in other words different evolutionary paths of the 
binary. Alternatively, it is possible that the diversity in 
T a arises from the diversity of the disk structure or the 
outburst mechanism. In order to investigate the origin of 
the T s diversity, detailed studies are required of the sub- 
groups which may be present in the shortest P or b regime, 
though the number of objects is small, apart from WZ Sge 
stars. 

In this paper we report detailed observations of a new 
DN having a short P orb . The object, OT J055717+683226 
(hereafter, J0557+68), was discovered on an image taken 
on 2006 Dec 16.6 (UT) by W. Klochr (Kloehr et al. 2006). 
Our time-series observations showed periodic modulations 
analogous to superhumps observed in SU UMa-type DNe 
(Uemura, Arai 2006). We also found that J0557+68 is a 
peculiar system in terms of its T s , which is much shorter 
than those of WZ Sge stars. 

As well as J0557+68, several objects in the shortest P orb 



regime have recently been discovered and well studied. 
The increasing sample size of short- P orb DNe now allows 
us to perform a quantitative classification of DNe. 

In the next section, we describe our observation equip- 
ment and image reduction. In § 3, we report the detailed 
outburst features of J0557+68. In § 4, we provide a new 
exploratory classification of DNe using a hierarchical clus- 
ter analysis. We also discuss the nature of J0557+68 and 
the origin of the diversity in T s in the shortest P olb regime. 
In the final section, we summarize our findings. 

2. Observations 

We performed optical and infrared observations of 
J0557+68 using TRISPEC, attached to the "KANATA" 
1.5-m telescope at Higashi-Hiroshima Observatory. 
TRISPEC is a simultaneous imager and spectrograph with 
polarimctry covering both optical and near-infrared wave- 
lengths (Watanabe et al. 2005). We used the imaging 
mode of TRISPEC with V, J, and K s filters. The effec- 
tive exposure times for each frame were 63, 60, and 54 s 
in the V, J, and K s bands, respectively. After making 
dark-subtracted and flat-fielded images, we measured the 
V magnitudes of J0557+68 using a comparison star lo- 
cated at R.A. = 05 h 57 m 29?105, Dec. = +68°29'18"88 
(V = 13.39). The J and K s magnitudes were measured 
using a comparison star at R.A. = 05 h 57 m 39?60, Dec. = 
+68°33'34"2 (J= 12.255 and K s = 11.937). We took the V 
magnitude of the comparison stars from the HIPPARCOS 
and TYCHO catalogues (Perryman, ESA 1997) and the 
J and K s magnitudes from the 2MASS catalog (Skrutskie 
et al. 2006). We checked the constancy of the comparison 
stars using neighboring stars and found that they exhib- 
ited no significant variations over 0.01, 0.02, and 0.05 mag 
in the V , J, and K s bands, respectively. Unfortunately, 
we could only obtain J-band images on JD 2454088 and 
2454089 during the first outburst due to a failure of a fil- 
ter wheel in the instrument. Heliocentric corrections were 
performed for each data set used for period analysis. 

We show example optical CCD images in figure 1. 
J0557+68 is marked with the black bars. The image in 
the left panel was taken on JD 2454088, during the out- 
burst. The right image was taken on JD 2454411, about 
11 months after the outburst when the object was at qui- 
escence. 

We also performed optical photometric observa- 
tions with other small telescope during the period of 
JD 2454087-2454114. The magnitudes with a clear fil- 
ter and without any filters were adjusted to the U-band 
system of TRISPEC/KANATA by adding constants. 1 
Table 1 details our observation log and instruments at 
each observatory. 

Using neighboring USNO B1.0 stars, we calculated 
the position of J0557+68 from 10 U-band images on 
JD 2454088 obtained with TRISPEC/KANATA. The 

The clear filter is made of the optical glass in a thickness com- 
parable to other filters. The total efficiency with the clear filter 
is slightly different from that without filters since the clear filter 
has its own transmission curve. 
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Table 1. Observation log and our equipment. 
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Site 


Telescope 


Camera 


Filter 


Date (+JD2454000) 


Higashi-Hiroshima 


1.5-m 


TRISPEC 


V, J, K s 


88, 89, 91, 92, 96, 97, 100, 105, 114, 


Astr. Obs. 








123, 125, 126, 127, 130, 131, 148, 151, 
152, 153, 155, 412, 543, 547 


Saitama 


25-cm 


SBIG ST-7XME 


clear 


87, 89, 92, 93, 97, 99, 543, 547 


Kyoto Univ. 


40-cm 


SBIG ST-9 


none 


101, 105, 109, 110, 111, 112, 113, 114 


Kagoshima Univ. 


1.0-m 




none 


101 




Fig. 1. V-band images of J0557+68 in outburst (left panel) 
and in quiescence (right panel) observed with the KANATA 
1.5-m telescope. The field of view is 5'x5'; north is up and 
east is to the left. The two black bars indicate J0557+68. 
The diffuse source indicated by the arrow in the right panel 
is a ghost image. 

average position is R.A. = 05 h 57 m 18?487, Dec. 
+68°32'27."01 with a systematic error of 0."18. Within 
the error of the position, a faint USNO B1.0 star with 
b = 19.04 and r = 19.23 could be observed, which is con- 
sidered to be the quiescent counterpart of J0557+68. A 
ROSAT X-ray source at 25" apart from J0557+68 could 
also be observed, which is probably the X-ray counterpart. 

3. Short-Period DN, J0557+68 

3.1. Overall behavior of the 2006 and 2008 outbursts 

The upper panel of figure 2 shows the light curve of 
an outburst of J0557+68 in 2006. The outburst con- 
tinued for 14 days from its discovery (JD 2454085.6; 
Kloehr et al. 2006), with an average fading rate of 0.10 ± 
0.01 magd^ 1 . These features are quite typical for super- 
outbursts in SU UMa-type DNe (Warner 1995b). This 
plateau phase was terminated by a rapid fading starting 
on JD 2454100. The object experienced a short rebright- 
ening on JD 2454109. The duration of the rcbrightening 
was ~ 5 days. Such single short rebrightcning is occasion- 
ally observed in short period SU UMa stars (e.g. Baba 
et al. 2000; Imada et al. 2006a) and WZ Sge stars (e.g. 
Ishioka et al. 2001b; Templeton et al. 2006). After the 
rebrightcning, the object again started a gradual fading. 
The object had almost returned to the quiescent state by 
JD 2454120. Our observation at quiescence (JD 2454412) 
yielded V = 18.89 ±0.04. The observed amplitude of the 
outburst is ~ 4 mag. We note that this is the lower- limit 
of the amplitude since we possibly overlooked the early 
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Fig. 2. Light curves of an outburst of J0557+68 in 2006 (up- 
per panel) and 2008 (lower panel). The abscissa and ordinate 
denote the time in JD and the magnitude, respectively. The 
filled circles are V-band and unfiltered CCD observations. 
The open circles and squares are J and A^-band observations, 
respectively. In the upper panel, the J and K B magnitudes 
are shifted by +0.5 and +1.0 mag, as indicated in the figure. 

phase of the outburst. 

We obtained simultaneous photometric data in the V 
and J-bands for two nights and in the V and .Kg-bands 
for four nights. These near-infrared data are shown in 
figure 2. Using these multicolor data, we show temporal 
color variations in figure 3. As shown in figure 3, both 
the colors remained ~ during the plateau phase, which 
is typical for DNe in outburst. In the rapid fading stage, 
a clear reddening of the object was observed in V — K s . 

The next outburst was observed in March 2008. We 
show the light curve of the 2008 outburst in the lower 
panel of figure 2. The object remained in a bright state for 
at least 4 days. The object was apparently fading at a rate 
of 0.07 magd -1 during the outburst. These characteristics 
indicate that the 2008 outburst is another superout burst. 
The time interval between the two superoutbursts is about 
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Fig. 3. Temporal evolution of colors. The abscissa and or- 
dinate denote JD and the colors of V — J and V — K B . The 
open and filled circles are V — J and V — K B , respectively. 
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Fig. 4. Example superhumps observed in J0557+68. The 
abscissa denotes the time in JD, in which integer parts are 
subtracted. The ordinate denotes the V-magnitude. Light 
curves on JD 2454087, 2454088, and 2454089 are shown. The 
light curve on JD 2454089 is shifted by +0.2 mag. 

480 d. We consider this to be a T s of J0557+68, which 
is much shorter than those of WZ Sge stars (~ lOyr) and 
rather comparable to those of normal SU UMa stars. 

3. 2. Superhumps 

During the 2006 outburst, we detected short-term pe- 
riodic modulations. Examples arc shown in figure 4. The 
modulations had saw-tooth profiles with decreasing ampli- 
tude with time. These features are typical characteristics 
for superhumps (Uemura, Arai 2006). Together with the 
features of the whole light curve, the detection of the su- 
perhumps established that J0557+68 is an SU UMa-type 
DN and the 2006 outburst was a superoutburst. 

We performed a period analysis for the superhumps 
using the phase dispersion minimization (PDM) method 
(Stellingwerf 1978). Before the PDM analysis, a linear 
fading trend was subtracted from the light curves dur- 
ing the superoutburst (JD 2454087-2454100). The light 
curves obtained with small telescopes have larger disper- 




16 18 20 22 
frequency (1/day) 

Fig. 5. Frequcncy-0 diagram defined by the PDM method 
(Stellingwerf 1978). The best candidate of the frequency is 
marked in the figure. 

sions than those obtained with the KANATA telescope. 
By binning these light curves into neighboring three-point 
bins, we obtained data with similar dispersions. We ex- 
cluded observations on JD 2454092 and 2454093 from the 
sample for our period analysis since the phase of modula- 
tions changed on these nights, as mentioned below. As a 
result, we obtained 620 photometric points. Using these 
data, we calculated defined in the PDM method for each 
frequency. The frequency-0 diagram is shown in figure 5. 
A strong signal can be seen at ~ 0.05341 d. 

Figure 6 presents the temporal evolution of the super- 
hump profiles. Short-term modulations observed after the 
superoutburst are also included in this figure. The light 
curves were folded by a period of 0.05341 d. As mentioned 
above, the supcrhump decreased in amplitude during the 
early phase of the superoutburst, as can be seen in fig- 
ure 6 (labeled as "JD2454087", "88", "89", and "91"). 
As the main hump weakened, the secondary hump ap- 
pears to become prominent. As a result, the hump had a 
double-peaked profile with the phase apparently inverted 
in several humps ( "92" , "93" , "97" ) . Just before and af- 
ter the rapid fading stage, the hump profiles became com- 
plicated, while hints of superhumps and their secondary 
humps can be seen ("99", "100", "101"). During the re- 
brightening, clear humps were observed, while their phase 
was shifted compared with the superhumps. They had si- 
nusoidal profiles rather than the typical saw-tooth profiles 
("109", "110"). 

To explore the temporal variation of the superhump pe- 
riod, we calculated O — C of the superhumps. The de- 
tailed procedure to determine the peak times of the su- 
perhumps is shown in Kato et al. (2009). We calculated 
O — C for them with a period of 0.05341 d and the epoch 
of HJD 2454087.2836, the first supcrhump maximum we 
observed. The O — C and the time of the supcrhump max- 
ima are listed in table 2. The O — C diagram is depicted 
in figure 7. In the table and figure, E denotes the cycle 
number of the superhumps. As described above, the phase 
of the prominent hump was inverted in the cycles 91, 92, 
and 94. O — C for these secondary humps are also listed 
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date 

+JD 2454000 



87 

[superoutburst) 



Table 2. O — C of superhumps. 



91 



92 



93 



97 



99 



100 

(rapid fading) 
101 



109 

(rebrightening) 
110 



-0.5 



0.5 
phase 



1.5 



E O — C Error Supcrhump peak 
(d) (d) (HJD-2400000) 






-0.0004 





0007 


54087 


2836 


1 


-0.0009 





0011 


54087 


3365 


15 


-0.0020 





0005 


54088 


0832 


16 


-0.0013 





0006 


54088 


1373 


17 


-0.0053 





0012 


54088 


1868 


19 


-0.0013 





0006 


54088 


2976 


20 


-0.0036 





0005 


54088 


3487 


34 


-0.0030 





0009 


54089 


0971 


35 


-0.0057 





0009 


54089 


1479 


36 


-0.0051 





0010 


54089 


2018 


37 


-0.0043 





0010 


54089 


2561 


38 


-0.0054 





0008 


54089 


3084 


39 


-0.0050 





0008 


54089 


3622 


49 


-0.0047 





0019 


54089 


8967 


50 


0.0007 





0012 


54089 


9555 


51 


-0.0027 





0011 


54090 


0055 


52 


-0.0041 





0042 


54090 


0575 


54 


-0.0014 





0031 


54090 


1670 


74 


0.0020 





0019 


54091 


2387 


91 


-0.0030 





0082 


54092 


1346 


92 


-0.0068 





0083 


54092 


1842 


94 


-0.0021 





0084 


54092 


2955 


(91) 


-0.0343 





0022 


54092 


1621 


(92) 


-0.0317 





0019 


54092 


2180 


(93) 


-0.0333 





0019 


54092 


2697 


(94) 


-0.0341 





0022 


54092 


3222 


109 


0.0085 





0048 


54093 


1147 


1 1 n 

11U 


n nno 

U.Ul LZ 


u 


UUOl 


o4Uyo 




143 


0.0041 





0021 


54094 


9264 


144 


0.0050 





0026 


54094 


9808 


180 


-0.0003 





0044 


54096 


8984 


183 


-0.0002 





0062 


54097 


0587 


185 


-0.0216 





0073 


54097 


1442 


257 


-0.0157 





0024 


54100 


9959 


258 


-0.0014 





0040 


54101 


0637 


259 


-0.0122 





0025 


54101 


1063 


260 


-0.0242 





0021 


54101 


1477 



Fig. 6. Phase-averaged light curves of superhumps. The ab- 
scissa and ordinate denote the supcrhump phase and the dif- 
ferential magnitude, respectively. The supcrhump phase was 
defined by an arbitrary epoch and the superhump period of 
0.05341 d. Temporal variations of the superhump profile arc 
shown from JD 2454087 (top) to JD 2454110 (bottom). Each 
light curve is shifted by 0.2 mag. 



in the table and figure. 

Figure 7 show a possible period change of superhumps. 
Kocn (2006) provides a method to test the significance of 
a period change seen in O — C diagrams. We applied this 
method to the O - C of J0557+68. Koen (2006) consid- 
ers two mechanisms which cause O — C variations, a ran- 
dom cycle-to-cycle jitter and a real period variation. Four 
models for O — C variations are then introduced: neither 
cycle-to-cycle jitter nor period change (Model 1), a signif- 
icant cycle-to-cycle jitter and no period change (Model 2), 
a significant period change and no cycle-to-cycle jitter 
(Model 3), and both a significant cycle-to-cycle jitter and 
a period change (Model 4). The goodness of fit of the 
models is evaluated with two information criteria (IC), 
the Akaike and Bayes IC (AIC and BIC). Probabilities 
for the models are finally calculated with each IC. We 
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Table 4. Best-fit parameters for the O — C model of 
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Fig. 7. O ~ C diagram of supcrhumps in J0557+68. The 
abscissa and ordinate denote the cycle number of supcrhumps 
and O — C in days, respectively. The dashed lines are the 
best-fit lines for O — C in the early and late stages (for detail, 
see the text). The open circles indicate the secondary humps 
prominent in cycles 91, 92, and 94. 



Table 3. Probabilities for the models in Koen (2006). 



Model* 


Pt 

J a 


Pt 




Early phase 


{0<E< 74) 


Ml 


0.02 


0.02 


M2 


0.20 


0.20 


M3 


0.64 


0.63 


M4 


0.14 


0.15 







Late phase (109 < E) 


Ml 


0.59 


0.41 


M2 


0.22 


0.28 


M3 


0.17 


0.22 


M4 


0.02 


0.09 



* For description of the models, see the text. 

t Probability calculated from Akaike information 
criterion. 

* Probability calculated from Bayesian information 
criterion. 



calculated the probabilities using the O — C observed in 
J0557+68 in an early phase of < E < 74 and a late phase 
of 109 < E. The results are listed in table 3. In this table, 
P a and P, denote the probabilities calculated from AIC 
and BIC, respectively. Model 3 and Model 1 have the 
highest probabilities among the four models in the early 
and late phases, respectively. 

The result of the Koen's test indicates that Psh of 
J0557+68 changed with time in an early phase, and then 
became constant in a late phase of the superout burst. The 



J0557+68. 


a 


2.87±0.56 x 10- b 


b 


-2.15 ±0.37 x 10" 5 


c 


1.78 ±5.28 x 10~ 4 


T 


108 ±7 


P 


-1.80 ±0.21 x 10" 5 


q 


3.02 ±0.70 x 10~ 2 


Plot 


10.8 ±2.2 x 10" b 


p f 


0.05324 ±0.00002 



* Period derivative in the early phase, 
t Period in the late phase. 



transition occurred in 74 < E < 109. Assuming a constant 
period derivative, these features can be modeled as fol- 
lows: 

(O - C% ar i y (P) = aE 2 + bE + c (0 < E < T), (1) 
(0-C) latc (E)=pE + q(E>T), (2) 
(0-C) early (T) = {0-C) late (T). (3) 

Equation (3) is needed in order to ensure the continuity 
of the O — C variation. The best-fit parameters and their 
uncertainties were calculated with a Bayesian approach 
in which probability densities of the parameters were esti- 
mated using a Markov-Chain Monte Carlo (MCMC) algo- 
rithm (Metropolis 1953; Gilks et al. 1996). Table 4 lists 
the obtained best-fit parameters for the observed O — C of 
J0557+68. We confirmed that the obtained T (= 108±7d) 
is consistent with that expected from the Koen's test 
(74 < E < 109). The model with those best-fitted pa- 
rameters is indicated by the dashed line in figure 7. The 
period derivative of the early phase (Pdot =Psh/Psh) and 
the period of the late phase (P2) are also included in the 
table. The positive period derivative means a period in- 
crease of supcrhumps in the early phase. 

According to Kato et al. (2009), SU UMa-type DNe gen- 
erally has three stages in term of Psh: stage A with a long 
Psh, stage B with a positive period derivative, and stage C 
with a short Psh- The O — C variation of J0557±68 can 
be interpreted as a transition from stage B to C. The du- 
ration of stage A is so short (1-2 d) that it was probably 
missed in the case of J0557±68. Kato et al. (2009) sug- 
gested that the minimum Psh, either P2 or Psh at the 
start of stage B, is regarded as a representative Psh of 
an object. The representative Psh of J0557±68 is, then, 
P 2 =0.05324 ±0.00002 d. 

To date, GW Lib has the shortest P or b among ordinary 
hydrogen-rich CVs except for V485 Cen and EI Psc, that 
is, P orb = 0.05332 ± 0.00002 d (Thorstensen et al. 2002b). 
VS 0329±1250 possibly has a shorter P or b since its super- 
hump period is quite short (P SH = 0.053394 ± 0.000007d; 
Shatter et al. 2007). As shown in table 4, J0557+68 is 
definitely one of SU UMa stars having the shortest P or b- 
V485 Cen and EI Psc are hydrogen-rich DNe having atyp- 
ically short P or b (~ 60 min). However, they are considered 
to have an evolutionary path different from those of CVs 
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Fig. 8. Upper panel: Phase-averaged superhump profiles in 
the V and J-bands on JD 2454088. The abscissa and or- 
dinate denote the superhump phase and magnitude, respec- 
tively. The open and filled circles are observations in the 
V and J-bands, respectively. Both the light curves are nor- 
malized to be readily compared. Lower panel: V — J color 
variation associated with the superhump. The dashed line 
indicates the color at the bottom of the super humps. 

having main-scqucncc secondaries (for details, sec § 1 and 
4). 

3.3. Simultaneous optical and near infrared observation 
of superhumps 

We successfully performed simultaneous V and J-band 
time-series observations on JD 2454088. Superhumps 
were observed in both bands. The phase-averaged profiles 
and colors are shown in figure 8. In the rising phase of 
the superhump, the 1^-band flux brightened more rapidly 
than the J-band flux. As a result, the bluest phase pre- 
ceded the superhump maximum. After the superhump 
maximum, the U-band flux decreased rapidly, while the 
decrease of the J-band flux was rather gradual. In addi- 
tion, we can see a secondary hump at a phase of ~ 0.7 
only in the J-band. The behavior in the J-band results 
in a rapid reddening, as can be seen in the lower panel 
of figure 8. The color during the secondary hump was 
redder than that at the bottom of the superhumps, as 
indicated by the dashed line in the figure. The reddest 
phase (~ 0.7) significantly preceded the bottom of the su- 
perhump (~ 0.0). 

It is noteworthy that the primary and secondary humps 
exhibited different behaviors in V — J. In the "blue" pri- 
mary hump, the observed feature can be explained by 
heating and subsequent cooling processes at the outermost 



part of the accretion disk. In superhumps, the heating 
mechanism is tidal dissipation (Osaki 1996). The prece- 
dence of the bluest phase suggests that the heating process 
finished before the superhump maximum. Then, the ob- 
ject reached the superhump maximum by an expansion 
of the low temperature region in the accretion disk. In 
the "red" secondary hump, the latter effect, namely, the 
expansion of a low temperature region, is probably more 
substantial. 

4. Discussion 

As shown in the last section, J0557+68 has a relatively 
short T s (480 d) and a quite short P or b ■ T s is much shorter 
than those of WZ Sge stars which are a dominant popu- 
lation in the shortest P or b regime. The characteristics 
of J0557+68 are reminiscent of V844 Her (Thorstensen 
et al. 2002b). As well as J0557+68 and V844 Her, recent 
studies have shown several objects in the shortest P or b 
regime having observational features distinct from those 
of WZ Sge stars, such as PU CMa (Kato et al. 2003b) and 
LL And (Kato 2004). On the other hand, there is no es- 
tablished classification of DNe in the shortest P or b regime 
that includes the above anomalies. In order to study the 
nature of ,10557+68, it is important to identify its pos- 
sible companions and search their common features. In 
this section, we present a new exploratory classification of 
DNe in the shortest P or b regime using hierarchical cluster 
analysis. We discuss not only the nature of J0557+68, 
but also the origin of the observed diversity in T s based 
on the identified clusters. 

4-1. Classification of DNe in the shortest P or b regime by 
hierarchical cluster analysis 

We list SU UMa-type DNe having short P orb (< 95min) 
in table 5. The table shows T s , outburst amplitude Am, 
and superhump period excess e = (Psh — Porb)/Porb- The 
table only includes confirmed SU UMa stars in which su- 
perhumps were observed during past superoutbursts. The 
objects and their parameters in the table were generally 
quoted from the Ritter & Kolb catalog ver. 7.9 (RKcat7.9; 
Ritter, Kolb 2003). Some of the parameters were updated 
by the references and recent superoutbursts shown in the 
table. The dates of the recent superoutbursts were re- 
ferred from the VSNET and ASAS-3 database (Kato et al. 
2004b; Pojmanski 2002). 

According to the disk instability model, T s and e can 
be indicators of the mass-transfer rate and the mass ratio 
of the binary, respectively (Osaki 1996; Patterson et al. 
2005a). Am has historically been considered to be a 
key parameter for dividing DNe into WZ Sge and normal 
SU UMa systems (Howell et al. 1995; Kato et al. 2001b). 

The parameters are plotted against P or b in figure 9. 
J0557+68 is indicated by the large open circle with its 
Psh- The anomalous feature of J0557+68 is evident: short 
T s and small Am, compared with the other short P or b sys- 
tems having 76min < P or b ^5 85min. According to the disk 
instability theory, the short T s suggests that the mass- 
transfer rate of J0557+68 is exceptionally high compared 
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with the systems having similar short P rb- 

Hierarchical cluster analysis was performed using the 
four parameters P or b, log(T s ), Am, and e. As the param- 
eters have different dimensions, for the cluster analysis the 
observed values were transformed to so-called "Z-scores" 
defined as (a — a)/a a , where a, a, and a a are the param- 
eter value, its average, and standard deviation, respec- 
tively. The sample consists of 34 objects selected from 
table 5 in which all four parameters are known. HO Del 
and AQ CMi were not included in the sample because 
their T s is uncertain. 

The calculation was performed with the pvclust 2 pack- 
age of R 3 . Using the pvclust function, we can estimate 
the confidence level (or probability value; p-value) for 
clusters via the multi-scale bootstrap re-sampling method 
(Shimodaira 2004). The bootstrap samples were gener- 
ated by randomly drawing N samples with replacement 
from the original sample, where N is the number of sam- 
ples. The p-value was approximated as the probability 
that the cluster is obtained in the bootstrap replicates of 
the dendrogram, shown in percentage. A high value indi- 
cates a high significance. In this paper, we only discuss 
clusters having p > 95 %. For the dendrogram classifica- 
tion method, we used Ward's method in which the within- 
class variance is minimized and the between-class variance 
is maximized (Murtagh, Heck 1987). This method is com- 
monly used for cases where the number of samples is small. 
It has the advantage of being able to derive small inde- 
pendent clusters, avoiding the need to form large chain- 
shaped clusters containing small clusters (Milligan 1980). 
Since our sample size is small, it is suitable for our anal- 
ysis. 

The obtained dendrogram is depicted in figure 10. We 
identified six clusters which have p-values higher than 95 
with a singleton cluster of DI UMa. The errors of the 
p- values were estimated to be less than 1 % in the range 
of p > 95 %. The seven clusters are labeled U SV", "DF , 
"EFT, tl WZ\ "V485", tl SU/WZ\ and "X" in the figure. 
The average parameters for each cluster are shown in fig- 
ure 11 as filled circles. We describe the characteristics of 
each cluster as follows: 

• Group SU consists of normal SU UMa stars. The 
minimum P or b is 86 min in this group. We consider 
that most SU UMa stars having P or b > 95 min also 
have the same nature as this group. 

• Group WZ includes so-called WZ Sge stars. This 
group is characterized by short P or b, long T s , large 
Am, and small e, as can be seen in figure 11. 

• Group SU/WZ is located between Group SU and 
WZ in all parameter spaces, as shown in fig- 
ure 11. All the systems in this group are listed as 
"large-amplitude SU UMa-type DNe" in Kato et al. 
(2001b), except for HV Vir. HV Vir has been con- 
sidered as a WZ Sge star since it exhibited early su- 
pcrhumps which are only seen in that group (Kato 
et al. 2001b). However, HV Vir is not classified into 

2 (http:/ /www. is. titech.ac.jp/~shimo/prog/pvclust/) 

3 (http://www.R-project.org) 
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Fig. 9. Supercycle, outburst amplitude, and superhump pe- 
riod excess of short period SU UMa-type DNe (P or b < 95 min) 
as a function of the orbital period. The open circles indicate 
systems categorized in a group having short P or b and short 
supercycles. The large open circles denote J0557+68 with 
Psh- The open triangles indicate V485 Cen and EI Psc. The 
filled circles indicate the other systems. 

Group WZ, but rather into SU/WZ in our cluster 
analysis because of its relatively short T s and large e 
compared with those in Group WZ. BC UMa is also 
a noteworthy member of this group. Maehara et al. 
(2007) reported the detection of early superhumps 
in this object, while BC UMa has an atypically long 
P or b for WZ Sge stars. They propose that BC UMa 
has an intermediate nature between normal SU UMa 
and WZ Sge stars. Group SU/WZ may be consid- 
ered as an intermediate evolutionary stage between 
Group SU and WZ. 
• Group V4-85 contains two peculiar systems, 
V485 Cen and EI Psc. They are indicated by open 
triangles in figure 9. As described in § 1, they have 
evolved secondaries and hence are proposed to be 
on a different evolutionary path from those of CVs 
having main-sequence secondaries (Augustcijn et al. 
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Fig. 10. Cluster dendrogram of short- P rb DNe. The vertical axis is a measure of the Euclidean dissimilarity. Shown are p-valucs 
for the clusters calculated with the multi-scale bootstrap re-sampling method (Shimodaira 2004). 



1996; Thorstensen et al. 2002a; Uemura et al. 2002). 

• Group ER is so-called "ER UMa stars", which are 
characterized by quite short T s (Kato, Kunjaya 
1995). In addition to ER UMa and VI 159 Ori, 
IX Dra is also considered to be in this group (Ishioka 
et al. 2001a). 

• Group DI is a singleton cluster of DI UMa. In gen- 
eral, this object is considered as a member of, so- 
called, "ER UMa stars" . Our cluster analysis di- 
vided ER UMa stars into Group ER and DI, prob- 
ably because of a very short P or b of DI UMa com- 
pared with those of ER UMa and VI 159 Ori. 

• Group X is a group established by our analysis for 
the first time. The members of this group are indi- 
cated by small open circles in figure 9. They have 
quite short P or b, although T s is about one order of 
magnitude shorter than those of " WZ" on average, 
as can be seen in figure 1 1 . Am is also much smaller 
than those of 11 WZ" and "SU/WZ*. e is significantly 
larger than those of " WZ' . 

Our classification method adequately reproduced the 
previously known sub-groups. In addition, we identified 
a new noteworthy group, Group X, which may be a key 
group for the study of CV evolution. 

4-2. On the nature of J 0557+ 68 and Group X 

J0557+68 presumably belongs to Group X because of 
the short T s and small Am. The result of the cluster 
analysis therefore predicts that J0557+68 has a relatively 
large e compared with those of WZ Sge stars. Using the 
average value of e in Group X (e = 2.2 ± 0.3), the P or b of 
J0557+68 is calculated to be 0.05209 ± 0.00015d. 

According to Podsiadlowski et al. (2003), CVs with 



evolved secondaries can have P m in shorter than that of 
CVs having main-sequence secondaries. Such objects have 
high mass-transfer rates before evolving to longer peri- 
ods compared with ordinary CVs, and finally evolve to 
AM CVn stars after passing P m i„. V485 Cen and EI Psc 
may be objects on this evolutionary channel (Uemura 
et al. 2002). As well as V485 Cen and EI Psc, objects 
on that evolutionary channel are expected to be found 
in the CV population with P or b > 76 min. Such progeni- 
tors of AM CVn and V485 Cen-like objects definitely have 
a comparable or higher mass-transfer rate compared with 
V485 Cen-like objects. As a result, they should have com- 
parable or shorter T s than those of ordinary objects. This 
might cause the observed diversity in T s in the shortest 
Porb regime (Thorstensen et al. 2002b). Possible candi- 
dates for such anomalies are ER UMa stars and objects 
in Group X. 

The mass-transfer rate in ER UMa stars should be 
one order higher than those of normal SU UMa stars in 
order to reproduce the quite short T s (Kato, Kunjaya 
1995; Osaki 1995a). If ER UMa stars are a progeni- 
tor of V485 Cen-like objects, the mass-transfer rate must 
rapidly decrease by one order of magnitude from P or b ~ 
90 min to ~ 60 min. The evolutionary path calculated by 
Podsiadlowski et al. (2003), however, shows a rather grad- 
ual decrease of the mass-transfer rate before P m i n . The 
small e of DI UMa, furthermore, indicates a small mass of 
the secondary. Thorstensen et al. (2002b) reported that 
the optical spectrum of DI UMa is typical for ordinary 
DNe at quiescence with none of the TiO absorption bands 
observed in V485 Cen and EI Psc. These observations 
are unfavorable for the scenario that ER UMa stars have 
evolved secondaries and are progenitors of V485 Cen-like 
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Fig. 11. As for figure 9, but for the averages of the parame- 
ters in sub-groups defined by our cluster analysis (for detail, 
see the text). The filled circles and error bars indicate the av- 
erages and standard deviations, respectively. The solid curves 
indicate the decision boundaries for each category, calculated 
by the support vector machine (SVM) method. 

objects. 

No evidence of evolved secondaries is seen either in the 
objects of Group X. The average e of the objects in Group 
X are smaller than those in Group V485, as can be seen 
in the bottom panel of figure 11. This implies that the 
objects in Group X have a secondary mass smaller than 
those of V485 Cen-like objects (Patterson et al. 2005a). 
Thorstensen et al. (2002b) reported that the optical spec- 
trum of V844 Her. a member of Group X, is also typical 
for ordinary DNe at quiescence. 

We noticed a peculiar outburst activity, observed in 
several objects in Group X. It is known that V844 Her 
experiences only a few normal outbursts and that most 
of the outbursts are superoutbursts (Kato, Uemura 2000; 
Oizumi et al. 2007). The disk instability model predicts 
that a system with a higher mass-transfer rate experiences 
more frequent normal outbursts (Osaki 1996). ER UMa 



stars actually show frequent normal outbursts, as pre- 
dicted. If V844 Her has an evolved secondary, the mass- 
transfer rate is expected to be high. However, the scarcity 
of normal outbursts appears to contradict disk instability 
theory. QZ Vir, a member of Group X, also exhibits few 
normal outbursts. Since QZ Vir is a bright source, it is 
less likely that a significant number of normal outbursts 
have been overlooked. In addition, no normal outbursts 
have been reported for J0557+68 based on the VSNET 
database (Kato et al. 2004b). 

The scarcity of normal outbursts in those systems is 
reminiscent of WZ Sge stars having a low mass-transfer 
rate. According to disk instability theory, the mass- 
transfer rate in WZ Sge stars is so low that the recurrence 
time of normal outbursts exceeds T s (Osaki 1995b). As a 
result, WZ Sge stars only show superoutbursts. It is pos- 
sible that the objects in Group X have low mass-transfer 
rates, like WZ Sge stars, while superoutbursts arc fre- 
quently triggered due to an atypically high viscosity in 
the quiescent accretion disk. In this case, the observed 
diversity in T s in the shortest P or b regime is not due to 
the diversity in the evolutionary path of objects, but due 
to the diversity in accretion-disk structures. 

The nature of Group X is still an open issue. Imada 
et al. (2006a) propose that the short -Porb 

DN SDSS 

J013701. 06— 091234.9 has a luminous and evolved sec- 
ondary since it shows significant TiO absorption and red 
infrared colors. As shown in table 5, the Am and e of 
SDSS J013701.06-091234.9 suggest that it is a possible 
member of Group X. This implies that a small portion of 
DNe having P orb > 76 min have evolved secondaries and 
hence are progenitors of AM CVn or V485 Cen-like ob- 
jects. Dynamical estimations of the secondary mass with 
radial velocity studies arc required for DNe in the shortest 
-Porb regime in order to explore the origin of the diversity 
in T s . Our cluster analysis provides high-priority objects 
for future studies, that is, members of Group X. 

4-3. Estimation of boundaries of groups with the support 
vector machine method 

The above seven groups of DNe were defined in four- 
parameter space. Now we estimate the boundaries of the 
groups in two-parameter spaces, that is, P rb and the 
other three parameters. This is useful for categorizing 
objects if only some of the parameters are known. 

The support vector machine (SVM) method provides 
a means to estimate the optimal decision boundaries of 
the groups (Vapnik 1998; Chen et al. 2005). The SVM 
constructs a linear classification between two groups by 
defining an optimal hyperplanc which separates members 
of the groups. The optimal plane is determined by maxi- 
mizing the margin between the opposing group members 
closest to the plane. The original SVM can be extended 
for non-linear classifications by using kernels to project 
the data into a higher dimensional space (Boser et al. 
1992). An improved SVM called the soft margin SVM 
can tolerate minor misclassifications (Cortes and Vapnik 
1995). In this paper, we used a soft-margin SVM with a 
radial-based kernel. 
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We estimated the boundaries of a group against the 
other groups using the result of the cluster analysis. 
The calculations were performed by the svm function (de- 
classification) in the el071 package of R. The width of an 
allowable margin around the separating plane is defined 
by a parameter "C". We determined the parameter "C" 
by maximizing the region of the groups and minimizing 
the overlapped areas between the regions. The calculated 
boundaries are indicated by solid curves in figure 11. As 
can be seen in figure 11, the boundaries fail to identify 
several groups: no boundary is given for Group ER and 
SU/WZ in the £-P or b plane, and the boundary of Group 
DI is given only for Am — P or b plane. This is because the 
number of the sample is small and the members of these 
groups are blended in the other group members in the 
s^Poib plane. 

In table 5, we show the group identification in paren- 
theses, for example "(SU)", for objects which were not 
used for the cluster analysis. These objects were classified 
following the boundaries estimated by SVM. In addition, 
objects were labeled "(WZ)" when i) early superhumps 
were detected and ii) the system experienced the rebright- 
ening phenomenon which characterizes WZ Sge-type su- 
peroutbursts (Kato et al. 2004a). Group SU, SU/WZ, and 
WZ and their possible members account for 70% and 95% 
of DNe with P or b < 86 min and 86 ruin < P or b < 95 min, 
respectively. Thus, we have confirmed that these three 
groups are the major population in those P or b regimes. 

5. Summary 

We observed J0557+68 during the outburst in 
December 2006. Our observation revealed that the object 
is an SU UMa-type DN having a quite short supcrhump 
period of P SH = 0.05324 ± 0.00002 d. The next superout- 
burst occurred in March 2008. The T s of this object is, 
hence, estimated to be 480 d, which is much shorter than 
those of WZ Sge stars. Using the cluster analysis, we 
found that a peculiar group characterized by short T s has 
P or b slightly longer than P m { n . J0557+68 probably be- 
longs to this group. While the nature of this group is still 
an open issue, its peculiar feature of a short T s is possi- 
bly due to an atypically high viscosity in quiescent disks 
or due to the evolutionary sequence being different from 
ordinary CVs having main-sequence secondaries. 

This work was partly supported by a Grand-in-Aid 
from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan (17684004, 17340054, 18740153, 
14079206, 18840032, 19740104). 
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Table 5. Parameters of DNe having a short period of P or b < 95 min. 
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Table 5. (Continued.) 
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